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We explore applications of the maximum entropy method (MEM) to de- 
termine properties of unstable particles using the four-dimensional 0(4) <f> 
\ theory as a laboratory. The spectral function of the correlation function of 

the unstable a particle is calculated with MEM, and shown to yield reliable 
■ results for both the mass of a and the energy of the two-pion state. Calcula- 

tions are also made for the case in which a is stable. Distinctive differences 
^vq \ in the volume dependence of the a mass and two-pion energy for the stable 

and unstable cases are analyzed in terms of perturbation theory. 

PACS number(s): 12.38.Gc, ll.15.Ha, 11.80.-m, 13.25.Jx 
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In contrast to many successes in understanding of stable hadron states by numerical 
simulations, unstable hadrons and hadronic decay processes are not understood well on the 
lattice. An important example in QCD is p — > tttt decay. There are a number of difficulties 
in dealing with such decay processes. One of the difficulties is computer power. In full QCD 
lattice simulations accessible today, the pion is so heavy that the p meson cannot decay. A 
more crucial problem is the difficulty first pointed out by Maiani and Testa [1]. Because the 
p meson and the tttt state in the isospin 1=1 channel have the same quantum numbers, the 
correlation function of the p meson behaves as a multi-exponential function corresponding 
to tttt states with various relative momenta as well as the p meson. In general, it is very 
difficult to decompose the states and extract the energy of each state from such a correlation 
function. 

In this article we explore application of the maximum entropy method (MEM) to extract 
the mass of unstable particle and the energies of states to which the unstable particle can 
decay. MEM has already been applied to extract the ground and first excited states of the 
stable mesons in QCD [2,3]. The masses of the two states were simultaneously extracted 
from the peak positions of the spectral function of the meson. This experience leads us to 
expect MEM to be useful for the multi-exponential correlation function of unstable particles 
such as the p meson. 

We employ the four- dimensional 0(4) 4 theory for the present exploratory study. It is 
an effective theory of QCD, containing the pion and a particle. The a can be made either 
stable or unstable by choice of parameters. In the non-linear formulation with the constraint 
<p a (x)(p a (x) = 1, the action is given by 



where <p a (x) = (tt 1 (x),(j(x)) for i = 1,2,3, and D(f>{x) = X^0(a; + A) + (fi( x ~ f 1 )- For 
m a > 2m w , the a particle is unstable and can decay to the 1 = two-pion state, similarly 
to the a meson in QCD. We examine the efficiency of MEM for extracting the energies of o 
and tttt states for both the unstable m a > 2m n and stable m a < 2m n cases. Very recently a 
similar study for the three-dimensional four-fermion model has been carried out by Allton 
et al. [4]. 

Numerical simulations are carried out at (/c, J) = (0.308,0.0012) and (0.30415,0.003), 
corresponding to m a jm^ m 3.7 and m a jm^ m 1.8, for several spatial lattice sizes in the range 
10 3 — 28 3 to investigate the volume dependence of the spectral functions and energies. The 
temporal lattice size is fixed to T = 64. Configurations are generated by the multi-cluster 
algorithm [5]. The periodic boundary condition is imposed in all directions. We perform 
0.6 x 10 6 - 1.2 x 10 6 iterations per simulation point (/c, J) and volume. 

We calculate the correlation function matrix [5,6] given by Cjj(r) = ((Oj(r) — 0«(r + 
l))Oj(0)) for i,j — a and tttt, where 0{ is the interpolating operator either for the o or 
the 1 = two-pion state with zero momentum, and the subtraction Oj(r) — Oj(r + 1) is 
made for eliminating the vacuum contribution. We apply MEM to the diagonal parts of the 
correlation function matrix Cu(t). 

In Fig. 1 we illustrate the correlation functions in the a and tttt channels for two volumes 
in the unstable case. The errors are estimated by the jackknife method with the bin size 
of 8000 to 16000 iterations. For the smaller volume the tttt correlation function exhibits a 
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change of slope at t ~ 10, from a larger slope corresponding to a to a smaller one of the mr 
state. The a correlator decays parallel to the latter, indicating dominance of the mr state 
for small volumes. For the larger volume, the trend is opposite, the mr correlation function 
showing a change of slope from the a particle to mi state at t ~ 20, while the mi correlator 
is dominated by the mi state. The decrease of "off-diagonal" contributions with increasing 
volume in the correlation functions can be understood from a perturbative analysis: the 
overlaps | (0 1 cr | 7T7t) | 2 and |(0|7r7r|cr)| 2 are proportional to l/[L 3 (m CT — E^) 2 ] where E nn is the 
two-pion energy. 

The spectral function f(uj) in the a and mi channels is defined in terms of the correlation 
function Cu(t) {i = a, nn) through 

L 3 Cii(T) = Jdufi(u)K(u,T), (2) 

where K{u,t) = e~™(l - e _aJ ) + e w(r ~ T) (l - e w ). In our MEM analysis the model function 
m{uj) is chosen as m(u) = A(u ) ■ Uq/lj 4 , where u is a reference point and A(u> ) is the 
perturbative a spectral function at a reference point u , which is given by 

AM = — /Zff . (3 ) 

■k[(uI -m 2 a + M r {uq)) + Mf(u )] K } 

The functions M R (u ) and Mj(u ) are defined by 

M R M = C[ ln M - ln (m a ) + 3 7<T (w )], (4) 
Mj(u; ) = Cni^M + 3/3 ff (wo)], (5) 



where 7 a (u;o) = PaM ]n[(l+p a (u ))/(l-p a (uo)] and p a (u ) = yjl - (4m 2 a /u$) 9(u -2m a ) 
for a = a and ir. The constant C in Eqs. (4) and (5) is given by C = 3Z n (rn 2 — m 2 ) 2 /32n 2 v 2 , 
where v = (a), and Z a (a = a, 7r) is the wave function renormalization factor for a and 7r, 
respectively. We apply the same model function to the a and hit correlation functions. 

We choose ojq = 2, and m a and Z a in Eq.(3) are fixed to the values estimated from inverse 
correlation functions in momentum space. The input parameters for the model function are 
compiled in Table I. The reconstruction is carried out in the region < u < 3. We choose 
Auj = 5 x 10~ 4 around the peaks corresponding to the a and tctt states to determine the 
energies accurately, and Auj = 10~ 2 in other regions. The number of data is taken as large 
as possible < r < T/2. We also check that the final results for the spectral functions do 
not depend much on the model function. 

In Fig. 2 we show the spectral functions for the a and tttt correlation functions, f a (ou) and 
f nn (ou), in the unstable case. Since energy eigenvalues are discrete on a finite volume, the 
spectral function is a sum of 5 functions. We indeed observe sharp peaks for both spectral 
functions, which can be identified as the mi state with zero momentum at the first peak and 
the a state at the second peak. The decrease of the peak height for the mr state in f a (uj) 
and for the a state in f nn (uj) agrees with the volume dependence of the correlation functions 
discussed above. In the figure for f a {uj) the position of peak expected for the mr state with 
momentum p = 2n/ L is indicated by a downward arrow. The absence of peak in our data 
implies that the overlap of a with this state is very small. 
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Fig. 3 shows the spectral functions in the stable case. For f nn (u) the a contribution 
decreases with volume similar to the unstable case. We observe only a single peak in f a (ou), 
indicating that | (0 1 1 7r7r ) | 2 is very small in this case. 

The a mass m a and the mr state energy E n7r obtained from the peak positions of the 
spectral functions as functions of spatial lattice size L are shown in Fig. 4(unstable case) and 
5 (stable case) by circles. For larger volumes m a from f nn (u}) and E nn from f a (u) suffer from 
large errors or are not available. This is because the overlaps | (0 1 cr 1 7r7r) | 2 and | (0 1 7r7r | ct) | 2 
decrease as the volume increases. 

We have seen in Figs. 2 and 3 that the a and txtx correlation functions are dominated 
by the a and nn states with zero momentum, and other states are negligible. In this 
case we can apply the diagonalization method [6] to the 2x2 correlation function matrix 
C(t) for extraction of the energy eigenvalues of these states. We diagonalize the matrix 
D(t, To) = C _1 / 2 (r ) C(r) C _1 / 2 (r ) at each r, where r is some reference time chosen to be 
t = in this work. The eigenvalue of D(t,t ) is given by \ u (t,t ) = K(W u ,t)/K(W u ,t ) 
with K{W V , t) = e~ rWv {l - e~ Wv ) + e w ^ T ~ T \l - e w »), where W v is the energy of the states 
v = a, txtx. 

The energies obtained by the diagonalization method are plotted by cross symbols in 
Figs. 4 and 5. The errors are estimated by the jackknife method with 8000 to 16000 elimina- 
tions. The results for m a and E wn are consistent with those with MEM, but the statistical 
error with the diagonalization method is much smaller. This is because the diagonalization 
makes full use of the 2x2 correlation function matrix while MEM utilizes on the diagonal 
element. 

For the application of the diagonalization method, it is essential that the dominant 
states in the correlation functions are known. Otherwise we need to consider the correlation 
function matrix of all possible states with the same quantum numbers, which is a difficult 
task. MEM is useful for identifying the dominant states, as exemplified with our example. 

Comparing Fig. 4 with Fig. 5 we find an essential difference in the volume dependence of 
the a mass and txtx energy between the unstable and stable cases. In the unstable case, the 
txtx energy increases and the a mass decreases as the volume increasing, while an opposite 
trend is seen in the stable case. 

The volume dependence of the a mass expected from perturbation theory is given by 

m a (L, m a , g R ) = m a + g R (Am a (L) - Am (T (oo)). (6) 
Here m a is the a mass for infinite volume, g R is the renormalized coupling constant, and 



4m a L 3 - a=7Vj(7 



+ A a C a (p) 



W a (p) 

with A a = 2, 6 for a = ir, a, D = 1 — 3(m 2 — m^/m 2 , and 
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Caip) W a (p)(ml-W*(p)y (8) 



where W a (p) = 2\Jm 2 a + 4 X)i=i sin 2 (j9j/2) with p { being the spatial momenta. 

We fit our results for the a mass obtained by the diagonalization method to Eq.(6), 
taking m a and g R as fit parameters and setting m w to the value obtained from the pion 
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correlation function at L — 28. The fit curves are plotted in Fig. 4 and Fig. 5 for each case, 
where the data at L — 10 are excluded from the fitting. The fit parameters and x 2 are 
compiled in Table II. The fits agrees quite well with the simulation results. We then realize 
that the different volume dependence between the unstable and stable cases originates from 
an opposite sign of C w (p) in Eq.(7) in the two cases. 

In order to understand the volume dependence of the 7T7T energy, we consider the scat- 
tering length do- It is related to the energy shift of the two-pion state through Liischer's 
formula [7] given by 



where c\ = —2.837297, c 2 = 6.375183. The results for a obtained from E nn calculated 
with the diagonalization method are tabulated in Table III in the column entitled "Simu. 
diago.". Here is fixed to the value at L — 28. The sign of a differ in the two cases, 
which reflects the difference in the volume dependence of E nn . 

We can estimate the scattering length in perturbation theory. From results for the 
perturbative phase shift [5] we obtain 



where R = ma/m^. To evaluate the right handside, we use the fit result for m a given in 
Table II. The coupling constant is obtained in two ways, either from the perturbative 
definition [5,8] denoted as g#(def.) = 3Z n (ml — m^)/v 2 or from the fit of the volume de- 
pendence of the a mass given in Table II denoted as ^(fit). For m w , Z n and v we use the 
results at L = 28 as before. 

In Table III the simulation results for ao are compared with the two estimates using 
the perturbative formula Eq. (10). While we cannot claim a precise agreement, we observe 
consistency in the value and sign obtained with simulation and perturbation theory for both 
the unstable and stable cases. Thus it is the factor of R 2 — 4 in the perturbative formula 
which leads to the the opposite volume dependence of the nn energy in the two cases found 
by the simulation. 

In this paper, we have investigated the efficiency of the maximum entropy method for 
study of unstable particle systems using the four-dimensional 0(4) 4 theory. We have 
demonstrated that the a mass and mi energy can be obtained from the a correlation func- 
tion alone. We have also explained that the difference in the volume dependences of the a 
mass and nn energy between the unstable and stable cases can be understood by pertur- 
bation theory. It is an advantage of MEM that only the single particle correlation function 
of the unstable particle is needed to analyze both the particle itself and the multi-particle 
decaying states. Furthermore it works even when the dominant states in the correlation 
functions are not known. We expect the MEM analysis to play a useful role in future stud- 
ies of unstable particles and decays in lattice QCD. 

This work is supported in part by Grants-in-Aid of the Ministry of Education 
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TABLES 



unstable case 

L m % m a Z a v 

10 0.1138(2) 0.961(1) 0.309(2) 0.956(1) 0.112921(7) 
12 0.1084(2) 0.9657(8) 0.337(2) 0.949(1) 0.127042(4) 
14 0.1053(2) 0.9667(6) 0.351(1) 0.944(1) 0.13351169(4) 
16 0.1041(3) 0.9674(6) 0.359(1) 0.938(1) 0.136818(1) 
18 0.1036(2) 0.9676(6) 0.362(1) 0.935(1) 0.138598(1) 
24 0.1024(2) 0.9678(8) 0.365(1) 0.928(1) 0.140646(2) 
28 0.1024(2) 0.9685(7) 0.365(1) 0.925(2) 0.1410983(6) 
stable case 

L m a Z CT v 

10 0.1945(2) 0.9734(4) 0.328(1) 0.967(1) 0.099507(1) 
18 0.1851(2) 0.9754(4) 0.3320(7) 0.951(1) 0.1094956(8) 
24 0.1844(2) 0.9735(5) 0.3344(6) 0.949(1) 0.1099262(1) 
28 0.1848(2) 0.9755(5) 0.3324(7) 0.945(1) 0.110067(2) 

TABLE I. The input parameters for the model function of the MEM in the unstable and stable 
cases. 



unstable stable 


m a 


0.3765(2) 0.3285(1) 


9r 


14(1) 9(1) 


X 2 /d.o.f 


. 0.23 2.5 



TABLE II. Fit parameters and x 2 /d.o.f. (degrees of freedom) of the volume dependence of the 
a mass. 





unstable 


stable 




a g R 


o 9r 


Simu. 


diago. 0.289(9) 


-2.49(19) 


Pert. 


g R (def.) 0.494(1) 19.1(2) 


-3.62(11) 17.8(1) 


Pert. 


fffl (fit) 0.361(3) 14(1) 


-1.98(23) 9(1) 



TABLE III. Scattering lengths for the the simulation (Simu.) and perturbative (Pert.) results. 
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FIG. 1. Correlation functions for the a and tttt in the unstable case. 
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FIG. 2. Spectral functions reconstructed from a (left line) and tttt (right line) correlation func- 
tions in the unstable case. 
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FIG. 3. Spectral functions reconstructed from a (left line) and ittt (right line) correlation func- 
tions in the stable case. 
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FIG. 4. Energies for a and tttt with the MEM and diagonalization in the unstable case. 
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FIG. 5. Energies for a and tttt with the MEM and diagonalization in the stable case. 
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